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ABSTRACT 

We report the Giant Metrewave Radio Telescope detection of HI 21cm absorption from 
the z ~ 3.39 damped Lyman-a absorber (DLA) towards PKS 0201+113, the highest red- 
shift at which 21cm absorption has been detected in a DLA. The absorption is spread over 
~ 115 km s- 1 and has two components, at z = 3.387144(17) and z = 3.386141(45). The 
stronger component has a redshift and velocity width in agreement with the tentative detection 
of Briggs et al. (1997), but a significantly lower optical depth. The core size and DLA cover- 
ing factor are estimated to be < 100 pc and / ~ 0.69, respectively, from a VLBA 328 MHz 
image. If one makes the conventional assumption that the HI column densities towards the 
optical and radio cores are the same, this optical depth corresponds to a spin temperature of 
T s ~ [(955 + 160) x (//0.69)] K. However, this assumption may not be correct, given that 
no metal-line absorption is seen at the redshift of the stronger 21cm component, indicating 
that this component does not arise along the line of sight to the optical QSO, and that there 
is structure in the 21cm absorbing gas on scales smaller than the size of the radio core. We 
model the 21cm absorbing gas as having a two-phase structure with cold dense gas randomly 
distributed within a diffuse envelope of warm gas. For such a model, our radio data indicate 
that, even if the optical QSO lies along a line-of-sight with a fortuitously high (~ 50%) cold 
gas fraction, the average cold gas fraction is low, (< 17%), when averaged over the the spatial 
extent of the radio core. Finally, the large mismatch between peak 21cm and optical redshifts 
and the complexity of both profiles makes it unlikely that the z ~ 3.39 DLA will be useful in 
tests of fundamental constant evolution. 
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1 INTRODUCTION 

While high redshift damped Lyman-o' absorbers (DLAs), the high- 
est HI column density systems seen in quasar spectra, are presum- 
ably the precursors of today's galaxies and the primary gas reser- 
voir for star formation, t heir size and struct ure have long been sub- 
jects of controversy (e.g. lWorfe et alj2005l) . Understanding the na- 
ture of a "typical" DLA as a function of redshift is one of the im- 
portant open problems in galaxy evolution. 

HI 21cm absorption studies of DLAs towards radio-loud 
quasars provide an estimate of the harmonic mean spin temperature 
of the absorbing gas, via a comparison between the 21cm optical 
depth and the HI column density obtained from the Lyman-a' line. 
This gives the distribution of gas in different temperature phases, 
which can be used to probe t he redshift evolution of ph ysical con- 
ditions in the absorbers (e.g. Kanekar & Briggs 2004). Further, a 
comparison between the redshifts of 21cm and optical metal-line 
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absorption (in a statistically large sample) can be used to study the 
evolu tion of a combination of fundamental constants 1 Wol fe et"al] 
1976). Extending 21cm absorption studies to the highest redshifts is 
interesting from both these perspectives. However, despite numer- 
ous searches, there are still only four confirm ed detections of 21cm 
absorption in z > 1 DLAs, out to z ~ 2.347 Jwolfe & Davislll979l : 
IWorfe&Briggslll98ll : IWolfeet all 1981 iKanekar et alj|2006l) . 

The z ~ 3.39 absorber towards PKS 0201+113 jWhite et al] 

Il993l) is one of the highest redshift DLAs known towards a 
radio-loud QSO and a lso has a very hig h HI column density 
(Nm ~ 1.8 x 10 21 cm~ 2 ; lElrison et al1 |200l|). Previous searches for 
21cm absorption have yielded conflicting results, with two tenta- 
tive detections (with s ubstantially different redshifts, line depths 
and velocity spreads; Ide Bruvn et al] 1 19961 ; iBriggs et al] 1997 ) 
and two non-detections at s imilar sensitivity dBriggs et al] 1997 ; 
Kanekar & Chengalur 1 19971) . all using different radio telescopes. 
High resolution Keck spectroscopy found no metal lines at either 
of the putative 21cm redshifts, suggesting that the 21cm detections 
might have been spurious tellison et alj200ll) . All the radio studies 
found evidence for a high spin temperature in the DLA, interpreted 
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Figure 1. GMRT HI 21cm spectra towards PKS 0201+113 in [A] 2001 July, [B] 2002 February, and [C] September 2003. at the original velocity resolution 
of ~ 7.2 km s _1 , with optical depth (—100 X t) plotted against heliocentric velocity, in km s -1 , relative to z = 3.38716. The main absorption component, close 

to zero velocity, is visible in all three spectra while the weaker component, at 70 km s -1 , can be seen in spectra [B] and [C], which have higher sensitivity. 

Panel [D] shows the final spectrum, after averaging the three spectra from [A]— [C]. The solid line in [D] shows the 2-Gaussian fit to the spectrum. The dashed 
vertical lines in each panel indicate the peaks of the two Gaussians of Table 1, while the error bar on the right side indicates the lcr noise on each spectrum. 



as a large fraction of warm HI. In contrast, the strong CII* lines in 
the Keck spectrum were used bv lWolfeetaUd2003l) to argue that 
the DLA contains a large fraction of cold HI. We report here on 
new deep Giant Metrewave Radio Telescope (GMRT) observations 
of the DLA in the redshifted HI 21cm line. 



2 OBSERVATIONS AND DATA ANALYSIS 

PKS 0201+113 was observed with the GMRT on multiple occa- 
sions, in 2001 January, 2001 July, 2002 February and 2003 Septem- 
ber. The GMRT FX correlator was used as the backend in all runs, 
with a 1 MHz bandwidth divided into 128 channels centred at the 
expected redshifted 21cm line frequency, giving a velocity resolu- 
tion of ~ 7.2 km s" 1 . The first three runs used the upper sideband 
of the correlator while the session in 2003 used the lower sideband. 
19 antennas were available for the observations in 2001 January, 24 
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Figure 2. GMRT Hanning-smoothed and re-sampled HI 21cm spectra towards PKS 0201+1 13 from [A] 2002 February, and [B] September 2003, with optical 
depth (-100 X r) plotted against heliocentric velocity, in km s _1 , relative to z = 3.38716. The spectra have a velocity resolution of ~ 14.4 km s -1 . The weaker 

component, at 70 km s , is more clearly visible in both panels. The dashed vertical lines in each panel indicate the peaks of the two Gaussians of Table 1. 

The error bar on the right side of each panel indicates the lcr noise on the spectrum. 



in July 2001, 28 in 2002 February and 27 in 2003 September. The 
varying UV-coverage is not a concern as PKS 0201+113 is unre- 
solved by the longest GMRT baselines. 

Observations of the standard calibrators 3C48 and 3C147 were 
used to calibrate the absolute flux density scale and the shape of 
the passband, while the compact source 0204+152 was used as the 
phase calibrator in all runs. The flux density of 0204+152 was mea- 
sured to be within a few percent of 6 Jy, the 327-MHz value quoted 
in the Very Large Array (VLA) calibrator manual, in all cases ex- 
cept in 2001 July, when it was found to be 5.62 ± 0.04 Jy. Note 
that the GMRT does not currently have online measurements of the 
system temperature; our experience indicates that the flux density 
scale is reliable to ~ 15% in this observing mode. 

Data from the different runs were converted to FITS and anal- 
ysed separately in classic AIPS, using standard procedures. After 
initial editing and gain and bandpass calibration, around 50 chan- 
nels (chosen to be away from the expected line location) were av- 
eraged together into a "channel-0" dataset and then used to obtain 
a continuum image of the field out to a radius of ~ 1 degree, well 
beyond the half power of the GMRT primary beam. Data from dif- 
ferent runs were imaged independently to test for source variabil- 
ity. The image was made by sub-dividing the field into 37 facets 
in all cases, to correct for the non-coplanarity of the array, and 
then used to self-calibrate the visibility data in an iterative man- 
ner until no further improvement in the image was obtained. This 
procedure typically converged after a few rounds of phase self- 
calibration, followed by two or three rounds of amplitude-and- 
phase self-calibration, with intermediate editing of corrupted data. 

The task UVSUB was next used to subtract all continuum 
emission from the multi-channel U-V data and the residual visi- 
bilities inspected for any radio frequency interference (RFI) on all 
baselines, using the task SPFLG to view the time-frequency plane. 
All detected RFI was edited out. The continuum-subtracted U-V 



data were shifted to the heliocentric frame (using the task CVEL) 
and imaged in all channels, and a spectrum then extracted from the 
image cube at the location of PKS 0201+113. Spectra were also 
extracted at the positions of other sources in the field to test for 
RFI. The data of 2001 January were found to be affected by low- 
level RFI and will hence not be discussed further; data from other 
epochs were found to be clean. The total on-source times were ~ 3, 
7 and 10 hours in 2001 July, 2002 February and 2003 September, 
respectively, after all editing. 

The 323 MHz flux density of PKS 0201+113 was measured 
to be 348 mjy, 417 mJy and 422 mJy in 2001 July, 2002 Febru- 
ary and 2003 September, from the GMRT images. The typical 
RMS noise in these images, away from bright sources, is ~ 0.35 - 
0.55 mJy/Bm. While the flux density in 2001 July is ~ 17% 
lower than in the other two runs, flux densities of other sources 
in the field (and the calibrator) were also lower by similar frac- 
tions; the GMRT data thus do not show any internal evidence for 
source variability. Ide Bruvn et al. Ul996h obtained a flux density of 
(350 + 10) mJy with the Westerbork Synthesis Radio Telescope 
(WSRT) in 1991 February, marginally consistent with the GMRT 
values (given the ~ 15% uncertainty in the GMRT flux scale); 
we were unable to compare flux densit ies of other sources in the 
0201+113 field for the WSRT dataset. iBriggs et all < fl997l) mea- 
sured (290 ±5) mJy in 1991 October with the VLA, significantly 
lower than the GMRT value. However, flux densities of other com- 
pact sources in the 0201+113 field in the VLA image are also 
lower (typically, by ~ 30%) than those in the GMRT image of 
2003 September. Given the consistency of the GMRT and VLA 
measurements relative to field sources, we feel that the differences 
in measured flux density are unlikely to be caused by variability 
in 0201+113 itself. The GMRT flux density scale is likely to be 
correct, as the flux density of the calibrator, 0204+152, is in good 
agreement with that in the VLA calibrator manual. We will hence- 
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Table 1. Parameters of the 2-Gaussian fit to the final GMRT spectrum. 
Notes : *Peak optical depth, using the GMRT flux density of 422 mJy. 
' Heliocentric frequency. ^Full width at half maximum. 



forth use a flux density of 422 mJy (from 2003 September) as the 
GMRT measurement. 

We also obtained a 328 MHz image of PKS 0201+113 with 
the Very Long Baseline Array (VLBA) in 2002, as part of a larger 
program to estimate the covering factors of quasars with foreground 
DLAs. The observations and data analysis will be described else- 
where; we merely note here that these data were analysed in the 
AIPS and DIFMAP packages, using standard procedures. The total 
on-source time was ~ 2 hours. 



3 SPECTRA AND RESULTS 

The final spectra obtained in 2001 July, 2002 February and 
2003 September are shown in panels [A]-[C] of Fig.[T] with optical 
depth (in units of - 100 x r) plotted as a function of heliocentric ve- 
locity, relative to z = 3.38716. The spectra are at the original veloc- 
ity resolution of ~ 7.2 km s -1 and have a root-mean-square (RMS) 
noise of 0.0057, 0.0043 and 0.0033 per channel, respectively, in 
units of optical depth (indicated by the error bars on the right side 

of each panel). The main absorption component, at 1.0 km s , 

can be clearly seen in all the spectra while a second, weaker, com- 
ponent is visible at 70 km s _1 in the higher sensitivity spectra of 

panels [B] and [C]. The latter two spectra are shown in Figs.[2]A] 
and [B], after Hanning-smoothing and re-sampling to a velocity 
resolution of ~ 14.4 km s , so that the weaker component can 
be seen more clearly. The doppler shift due to the Earth's motion 
between the three runs (e.g. ~ 42 km s _1 or ~ 45 kHz between 
2002 February and 2003 September) is far larger than the spectral 
resolution; both features are hence likely to be real. 

Fig. [TjD] shows the final GMRT spectrum towards 
PKS 0201+1 13, obtained from a weighted average of the three op- 
tical depth spectra in panels [A]-[C] of Fig. Q] The RMS noise 
on this spectrum is 0.0024 per 7.8 kHz channel (in optical depth 
units) while the equivalent width of the absorption is J T j s dV = 
(0.714 ± 0.017) km s -1 . The solid line in Fig.[T]T>] shows a two- 
Gaussian fit to the spectrum; the parameters of the fit are sum- 
marised in Table 1 . 

The 328 MHz VLBA image of PKS 0201+113, shown in 
Fig. [U has an angular resolution of ~ 78 x 38 mas 2 (i.e. ~ 588 x 
290 pc 2 at z ~ 3.39, the DLA redshift). The source is well-fit by a 
single elliptical Gaussian component; we measure a core flux den- 
sity of ~ 292 mJy and a deconvolved source size of 17.6 x 6.6 mas 2 
(i.e. ~ 132 x 50 pc 2 at z ~ 3.350). A circular Gaussian model 
gives the same core flux density and an almost identical x~ f° r the 
fit; the deconvolved source size is then ~ 11.6 x 11.6 mas 2 (i.e. 



1 We use the standard ACDM cosmology throughout this paper, with fl„, 
0.27, H A = 0.73 and H = 71 km s -1 Mpc~' 
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Figure 3. VLBA 328 MHz continuum image of PKS 0201+113, with a 
resolution of 78 X 38 mas. The contours are at 7.5 X (-1, 1,2, ...,64) mJy. 



~ 87.5 x 87.5 pc at the absorber redshift). We note that attempts 
to fit an additional component to the marginal north-west extension 
did not significantly alter the core flux density. 



4 DISCUSSION 

In this section, results of the present GMRT observations are ini- 
tially compared with earlier 21cm studies of the z ~ 3.39 DLA. 
The gas kinematics of the DLA are briefly discussed before we 
compare the 21cm and optical low-ionization metal profiles, as a 
probe of fundamental constant evolution. We next derive the spin 
temperature of the absorbing gas under the standard assumption 
that the average HI column density towards the radio core is the 
same as that measured towards the optical QSO from the damped 
Lyman-d' line. The comparison between metal line and 21cm pro- 
files finds evidence for small-scale structure in the absorbing gas 
across the radio core, suggesting that the average column density 
against the radio core need not be the same as that towards the opti- 
cal QSO. Finally, we discuss the HI temperature distribution in the 
light of the above results and two-phase models, and show that the 
gas must be predominantly warm in the absorber. 



4.1 A comparison with earlier 21cm absorption studies 

As noted in the introduction, the z ~ 3.39 DLA has been the 
focus of a number of 21cm absorption stud i es in the pas t , with 
very different results |de Bruvn et"ai] 1 1 9961 : iBriggs et alj 1 19971 ; 
Kanekar & Chengalur 1 1997 ). Details of the different observations 
are listed in Table 2, with the line depths quoted in both optical 
depth and flux density, while Fig.|4]shows a visual comparison be- 
tween the GMRT spectrum and the earlier results. Note that the 
comparison in the figure is in flux density units (rather than in opti- 
cal depth), as the continuum flux density estimates at Arecibo and 
the ORT are unreliable due to the very large beamwidths and the 
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WSRT 


1991 


350 


9.0 


7.0 


8.5 ±2.0 


30.0 ± 7.0 


3.38699(3) 


9±2 
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Arecibo 


1993 


318 


4.5 


4.1 


3.7 ± 0.8 


11.6 ±2.0 


3.38716(7) 


23 ±5 
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VLA 


1991 


290 


11.1 


3.5 


1.5 ± 1.2 


4.2 ±3.4 
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ORT 


1997 




5.6 


3.7 


1.3 ±0.6 


4.7 ± 2.2 






3 


GMRT 


2001 - 2003 


422 


7.4 


1.0 


1.76 ±0.20 
0.88 ±0.15 


7.43 ± 0.84 
3.71 ±0.63 


3.387144(17) 
3.386141(45) 


18.6 ±2.6 
35.1 ±7.4 


4 



Table 2. A comparison between the 21cm spectra from different telescopes. 'The quoted line and optical depths are from the Gaussian fits to the lines for 
the WSRT, Arecibo and GMRT spectra. For the VLA and ORT non-detections, the quoted line and optical depths are the values at 323.764 MHz (ORT) and 
323.770 MHz (VLA), after smoothing the spectra to resolutions of ~ 16.8 km s _I (ORT) and ~ 22.2 km s -1 (VLA), respectively, and resampling. The ORT 
optical depth uses the WSRT flux density of 350 mjy, as the ORT observation did not have a reliable simultaneous flux density measurement. Note that the 
WSRT and VLA values have n ot b een scaled to the GM RT flux density value of 422 mjy but are on the original flux scales. References : (l) lde Bruvn et alj 
GUI); (2) lBriggsetai]ll997l) ; (3) lKanekar & ChengaluJ dl997ri ; (4) this paper. 
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Figure 4. A comparison between the final GMRT spectrum (solid line) 
and earlier Arecibo, WSRT, VLA and ORT results, with flux density (in 
mjy) plotted against heliocentric velocity, in km/s, relative to z. = 3.38716. 
The dashed and dotted lines show the Gaussian fits to the earlier tentative 
Arecib o and WSRT detections of iBriggs et alj d 19971) and Ide Bruvn et al J 
d 19961) . respectively. The solid and open circles (and lcr err ors) indicate 
the values at 323.764 MHz (ORT) and 323.770 MHz (VLA) dBriggs et alj 
ll997l ; lKanekar &"c hengalur 1997) respectively, after smoothing to resolu- 
tions of ~ 16.8 km s" 1 (ORT) and ~ 22.2 km s _1 (VLA). The VLA and 
WSRT line depths have been scaled to the GMRT flux density of 422 mjy 
for the comparison. The 1<t errors on the Arecibo, WSRT and GMRT fits 
are shown by the error bars to the right of each telescope name. See text for 
discussion. 



associated confusion at the low observing frequency. The VLA and 
WSRT line depths have been scaled to the GMRT flux density of 
422 mjy for the comparison. 

It is clear from Fig. 0] and Table 2 that the VLA and ORT 
non-detections are consistent with the GMRT detection within the 
lcr errors; note that both spectra showed weak absorption (within 
the noise) close to the frequency of the main GMRT component 
(~ 323 .76 MHz). Next, the fit to the Arecibo feature of lBriggs etail 
dl997h is in excellent agreement with the stronger GMRT 21cm 
component in both redshift and FWHM. However, the peak opti- 
cal depth is lower in the GMRT spectrum, discrepant at ~ 2.4cr 
level. Similarly, the Areci bo peak optical depth diffe rs at ~ 2tr 
level from the ORT limit dKanekar & Chengalur|[l997h and is only 
in ma rginal agreement with the VLA non-detection Briggs et"al1 
1 1997b . A plausible cause of the optical depth discrepancy is an in- 
correct continuum flux density in the Arecibo data, a not-unusual 
problem with single-dish measurements. However, the line depth, 
in flux density units, in the Arecibo spectrum is also discrepant 
from the GMRT measurement at ~ 2cr level, implying that the dif- 
ference in optical depths is not solely due to an err or in the contin- 
uum flux density. As noted bv lBriggs et al. I d 19971) . low-level RFI 
was present throughout the Arecibo spectrum (again, not uncom- 
mon for single-dish observations), resulting in non-Gaussian be- 
haviour and possible systematic errors. This could well result in an 
incorrect determination of the Arecibo line depth and, hence, an 
over-estimate in the optical depth. 

The only signific antly discrepan t result in Fig.|4]and Table 2 is 
the WSRT feature of Ide Bruvn et al.l d 19961) . which is inconsistent 
with all the other data. The WSRT feature is unre solved and was 
detec ted at only ~ 4.5cr level in one polarization dde Bruvn et al.l 
1996); further these data have the lowest sensitivity of all the mea- 
surements. Given the difficulty of these low frequency observations 
in the presence of RFI, the relatively low statistical significance and 
the disagreement with all the other observations, the WSRT feature 
appears likely to be an artefact. 

One possibility is that the above differences between the 
WSRT, Arecibo and GMRT features might be due to real changes 
in the optical depth of the 21cm absorption, as has been seen 
in two low z DLAs, towards AO 0235+164 dWolfe et al .11 19821) 
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and PKS 1127-145 dKanekar & Chengalunl2001ah . However, the 
GMRT spectra show no internal evidence for optical depth vari- 
ability between the different observing epochs (from 2001 to 2003) 
and are also consistent with the VLA and ORT spectra (of 1991 and 
1997), making optical depth variability an u nlikely possibil i ty. We 
hence conclude that (1) the WSRT feature of de Bruvn et al. I fl996h 
is likely to be an artefact, (2) the GMRT, ORT, VLA and Arecibo 
spectra are consistent, if one assumes an incorrect estimate of the 
line depth in the Arecibo spectrum (note that the Arecibo spec- 
trum is only discrepant at ~ 2a level even without this assumption, 
with the FWHM and redshift in good agreement), and (3) there is 
no strong evidence for variability in the 21cm optical depth in the 
z ~ 3.39 DLA. Of course, it is obvious from the RMS values of 
Table 2 that the second, weaker, feature seen in the GMRT spec- 
trum is far too weak to have been detected in any of the earlier 
observations. 



4.2 Gas kinematics 

The 21cm absorption profile of Fig.[TjD] has a total velocity spread 
(between nulls) of ~ 115 km s" 1 , the largest 21cm velocity width 
seen in any z > 2 DLA. In fact, this width is comparable to values 
seen in low z DLAs identified as bright spiral gala xies (for compar- 
ison, see Table 3 in iKanekar & Chengaludl2003l) . The only low z 
DLA with a 21cm velocity width larger than ~ 100 km s _I and not 
associated with a spiral disk galaxy is the z ~ 0.3127 absorber 
toward s PKS 1127-145 dLane et all 1 1991 IKanekar & Chengaluil 
l2001ah . Interestingly enough, Ellison et al. ( 2001 ) find the gas kine- 
matics of the z ~ 3.39 DLA towards PKS 0201+113 to be chaotic, 
with no evidence for an edge-le ading asymmetry in the op tical lines 
(a signature of a rotating disk; IProchaska & Wolfein"997l) . Instead, 
the low-ionization metal lines show a large number of spectral com- 
ponents over a velocity spread of ~ 270 km s _1 . Similarly, both 
21cm components in the z ~ 3.39 DLA are quite wide (see FWHMs 
in Table 1), perhaps due to multiple blended sub-components. 



4.3 Fundamental constant evolution 

Comparisons between the redshifts of different spectral transitions 
can be used to probe evolution in the values of the fundamental 
constants such as the fine structure constant a, the electron-proton 
mass ratio /j = m e /m p , the proton gyromagnetic ratio g„, etc. A 
variety of approaches employing different transitions have been 
used , constraining changes in the above constants to be < 10~ 5 
(e.g.lMurphv et alj2003l ; isrianand et alj2004lKanekar et al]2005t 
Rein hold et alj|2006t) . One such technique uses optical metal lines 
and the 21cm line to probe changes in the quantity X = g p fia 2 
dWolfeetalJll976l) . However, the velocity structure can be very 
different in the optical and radio lines, making it difficult to iden- 
tify a single absorption redshift for the comparison, especially in 
cases where the profiles are complex, with multiple spectral com- 
ponents. While one might compare the redshifts o f the strongest ab- 
sorpti on component in the 21cm and metal lines l lTzanavaris et al.l 
2006), there is no physical reason why the strongest optical and 
21cm absorption should arise in the same spectral component; this 
depends on local physical conditions in the absorbing clouds. The 
"strongest-component" approach thus introduces an extra source of 
systematic error that could lead to a bias in favour of a detection in 
small absorption samples. Conversely, while a comparison between 
the redshifts of "nearest-neighb our" absorpt ion components in the 
21cm and optical lines ( IKanekar et alj [2006) might be viewed as a 



test of the null hypothesis, it has the opposite bias, i.e. towards a 
non-detection of evolution. 

The z ~ 3.39 DLA towards PKS 0201+113 is the highest 
redshift absorber in which 21cm and optical lines have both been 
detected and is thus an obvious candidate for studies of changes 
in X = gp/ua 2 . Unfortunately, both the 21cm and optical profiles 
are complex, extending over > 100 km s _I , with multiple compo- 
nents. Fig. [5] shows the GMRT 21cm profile (histogram) overlaid 
on the unsaturated FeII-/ll 122 and CIP - /H33 5 lines detected in the 
Keck-HIRES spectrum o f lEllison et all d200lh . It can immediately 
be seen that even the strongest optical component is not uniquely 

determined; the absorption at 65 km s -1 is stronger in the Fell- 

A.1 122 transition while that at 25 km s is stronger in the CII* 

line. We will use the CII* transition f or the comparison, as this is 
expected to primarily arise in cold gas (Wolf e etal .1120031 ; hereafter 
W03), which also gives rise to the 21cm absorption. However, simi- 
lar results are obtained on using other low-ionization metal species. 

The CII* absorption has three components, at za - 
3.386821(7), z B = 3.386297(10) and z c = 3.386541(10), in or- 
der of decreasing strength (the redshifts are indicated by the arrows 
in Fig.[5]Fi]), while the 21cm profile has two components, the main 
one at z, = 3.387144(17) and the secondary at z 2 = 3.386141(45). 
A strongest-component comparison between the 21cm and CII* 
lines (i.e. between za and zi) gives [AX/X] = [Az/(1 + z)] = 
(-7.4 ± 0.4) x 10~ 5 , where Az = Zai, - z 2 i and z = fei + Zcn.)/2. 
If this value of [AX/X] is due to changes in [Aa/a] or [Afi/jj.], 
it would require [Ace/a'] or [A/i/fi] to be nearly an order of mag- 
nitude larger than value s obtained from other methods probing a 
similar redshif t rangt El dMurphv et alj|2003l; ISrianand et al1l2004l; 
iReinhold etal] l2006h . Further. iTzanavaris et alj d2006l) obtained 
[AX/X] = (0.63 + 0.99) x lO" 5 over 0.23 < z < 2.35 from a 
strongest-component comparison between 21cm and optical lines 
in nine absorbers, again significantly smaller than the present value. 
This suggests that the difference between the redshifts of the 
strongest components towards PKS 0201+113 stems from intrin- 
sic differences between the radio and optical sightlines, discussed 
further in the next section. 

Alternatively, one might compare the redshifts of the weaker 
21cm feature and its nearest-neighbour CII* component (i.e. com- 
ponent B); formally, this yields [AX/X] = (+3.6 ± 1.0) x lO" 5 . 
However, these CII* and 21cm components have very different ve- 
locity widths (FWHMs of ~ 13 and ~ 21 km s" 1 , respectively), 
suggesting that at least part of the absorption does not arise in the 
same gas. A nearest-neighbour comparison between these optical 
and 21cm components is thus unlikely to yield reliable results. We 
hence conclude that, despite its high redshift, the z ~ 3.39 DLA is 
unlikely to be useful for the purpose of probing fundamental con- 
stant evolution. 

4.4 Physical conditions in the absorber 

The HI 21cm absorption in DLAs has generally been regarded 
as arising in a two-phase medium, with both cold (CNM) and 
warm (WNM) gas, similar to that seen in the Galaxy. Observa- 
tional support for such a model has been found in th e case of two 
low z DLAs dLaneetal.ll2OO0l; IKanekar etal] 1200 lh . In the con- 
text of such two-phase models, t he fact that high z DLAs generally 
have high spin temperatures (e.g. lde Bruvn et alTl 996; Caril li et al.l 

2 Note that chang es in g D are expected to be much sm aller than those in a 
and \i ICalmet & Fri tzsch 2002; lFlambaum et al.l2004l) 
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Figure 5. A comparison between the GMRT HI 21cm profile (histogram) and the [A] Fell Al 122 and [B] CII* /11335 Keck-Hires profiles (solid squares). The 
spectra are in arbitrary units and plotted against heliocentric velocity, in km s , relative to z = 3.38716. The dashed vertical lines in each panel indicate the 
redshifts of the two 21cm components from the 2-Gaussian fit. The three vertical arrows in [B] indicate the redshifts of the three components (A, B and C) 
obtained from the fit to the CII* line. See text for discussion. 



1996; Kan ekar & Chengalurl 1200 lbl l2003h implies that the bulk 
of the gas is in the warm phase. Recently, however, W03 used 
the observed CII* absorption in DLAs as a tracer of the cool- 
ing rate to argue that DLAs with strong CII* lines must have a 
substantial fraction o f their gas in the CNM phase [see, however, 
ISrianand et alj d2005l) l. W03 ruled out models in which all of the 
observed CII* absorption arises in the WNM, finding that strong 
CII* absor ption requi r es tha t the line-of-sight includes some CNM 
[see also lWolfe et al.1 (2004)]. They further argued that the multi- 
ple lines of evidence sugge sting that HI in high z DLAs consists 
primarily of warm gas (e.g. lLisztll2002l ; iNorman & Spaans|[T997l ; 
iKanekar & Chengalurl I2OO3I) are inconclusive, and suggested in- 
stead that a two-phase model with a CNM fraction of 50% towards 
the optical QSO is consistent with the observed bolometric lumi- 
nosity in DLAs. 

As noted by W03, the z ~ 3.39 absorber towards 
PKS 0201+113 is one of only two high z DLAs that has both 
strong CII* absorption (suggesting a high CNM fraction in the W03 
model), as well as a high 21cm spin temperature. W03 resolved 
this apparent contradiction by arguing that the background radio 
core must subtend a linear diameter of > 40 pc at the DLA, while 
the CNM "clouds" in the absorber have linear sizes of < 10 pc. 
They suggested that the strong CII* absorption could be consistent 
with a lack of 21cm absorption if the CII* absorption comes from 
a number of small CNM clouds lined up against the AU-sized op- 
tical QSO. We revisit this issue in the light of the new observations 
presented in the preceding sections. 

A convenient starting point for this discussion is the general 
relation between 21cm optical depth m, total HI column density 
A'hi (in citT 2 ) alon g the path, and spin temp erature T s (in K) in the 
optically thin limit dRohlfs & Wilsoij2004l) . 



1.823 x 10 18 T. 



where the integral is over velocity, in km s -1 , The spin tempera- 
ture thus obtained is the column-density-weighted harmonic mean 
of the T s values of different HI phases along the path. As such, it 
does not give the kinetic temperature of any individual phase but 
instead contains information on the distribution of gas in different 
temperature phases along the line of sight. 

Next, the total flux density of the source measured by a short- 
baseline interferometer such as the GMRT or the VLA may be 
much larger than the flu x density that is covere d by the foreground 
absorbing "cloud" (e.g. lBriggs&Wolfelll983l) . This is especially 
true at the high redshift of the present observations; for example, 
the GMRT synthesized beam in the dataset of September 2003 is 
~ 12.2" x 10.0", subtending spatial dimensions of ~ 100 x 75 kpc 
at z ~ 3.39, far larger than the size of a typical galaxy. In such 
cases, the measured 21cm optical depth T obs will be smaller than 
the "true" 21cm optical depth, t 2 i. If only a fraction / of the radio 
continuum is covered by the foreground absorbing gas, the "true" 
21cm optical depth is given by r 2 i = (1//) * T obs . For a damped ab- 
sorber, where the HI column density is measured from the Lyman-a 
line towards the optical QSO, we then have 



N m = 1.823 xl0 15 (T s //) T o6l dV 



(2) 



r„dV 



(1) 



/ is referred to as the 21cm covering factor of the DLA. 

In addition to the assumption that the 21cm absorption is opti- 
cally thin, Eqn. [2] contains two critical unknowns, the DLA cov- 
ering factor /, and whether the HI properties along the line of 
sight to the optical QSO can be applied to the much larger line 
of sight towards the more-extended radio quasar core. The cover- 
ing facto r has been the subject of much discussion in the litera- 
ture (e.g. iBriggs & Wolfdfl983l ; ICurran et al.ll2005l) . In principle, 
the "true" 21cm optical depth can be directly determi ned by us- 
ing V LB I techniques to map the 21cm absorption (e.g. lLane et al.l 
2000); however, this is usually impossible due to the fairly poor 
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frequency coverage and sensitivity of VLBI receivers. The alterna- 
tive approach is to use VLBI observations at nearby frequencies to 
measure the fraction of flux density arising from the compact ra- 
dio core. Even such measurements are rarely carried out in practice 
due to the tec hnical difficu lties in such VLBI studies at low fre- 
quencies (e.g. lBriggj|l983h . On the other hand, it is difficult, even 
in principle, to determine whether the HI distribution is uniform 
across the exte nded background radio source [see, for example, the 
discussions in Ide Bruvn et alj d 19961) and W03]. Given that these 
two assumptions are generally untested in high z DLAs, the case 
of PKS0201+113 is doubly unusual in that we have meaningful 
VLBI constraints on the size of the radio core, and, as discussed 
below, also find observational evidence for small-scale structure in 
the HI 21cm absorbing gas. 

For ease of comparison with earlier work, however, we first 
estimate the spin temperature of the z ~ 3.39 DLA with the usual 
assumption that the HI gas is uniformly distributed across the ra- 
dio core. If the 21cm absorber is a uniform slab covering the entire 
VLBI core, the covering factor / is S VL ba/S gmrt ~ (292/422) ~ 
0.69. Further, N m = (1.8 ± 0.3) x 10 21 cm" 2 dEllison et alj|200ll) . 
and J T ohs dY = (0.714 ± 0.017) km s _1 , which gives T s ~ [(955 ± 
160) x (//0.69)] K. This is in the "high" range of T s values in 
DLAs 1 Kaneka T&Chengalurl jpoj) . We note that the above T s 
value is significantly l ower than the estimate T s > 3290 K in 
iKanekar & Chengalurl J2003I) . This is because the latter assumed 
/ ~ 1 and a thermal velocity spread (FWHM = 20 km s _I ) which 
is much smaller than the true velocity spread of ~ 115 km s . 

We next consider the assumption that the gas is uniformly 
distributed across the radio core. If the radio source were indeed 
very compact (< 10 pc), strong interstellar scintillation would be 
expected at the low GMRT observing frequencies, giving rise to 
a variable source flux density. However, the GMRT observations 
found no evidence for flux density variability over time-scales of a 
few years. It thus appears likely that the radio core is at least a few 
tens of parsecs in size at the redshifted 21cm frequency, far larger 
than the size of the optical core. Our VLBA observations cannot 
rule out this possibility as they only require the radio core to be 
smaller than ~ 100 pc, the size of the fitted circular or elliptical 
Gaussian models. 

The large difference in the sizes of the optical and radio cores 
implies that a comparison between the profiles of low ionization 
metal lines and the 21cm line would reveal the presence of small- 
scale velocity structure in the absorber over scales of tens of par- 
secs. This is because the low ionization metal lines trace the kine- 
matics of the gas along the line of sight to the AU-sized optical 
QSO, while the 21cm line traces the kinematics along the much- 
larger (> 40 pc) region lying in front of the radio core. No metal- 
line absorption can be seen in Figs.[5]A] and [B] at z ~ 3.387144, 
the redshift of the stronger 21cm component. As noted in Sec- 
tion l4.3l it is not possible to account for the large observed redshift 
difference between the 21cm and optical l ines by fundamental con- 
stant evolution. Further, as pointed out by Ellison et al. (200lth the 
21cm peak redshift is at the edge of the absorption seen in even the 
saturated OI- /11302 and CII-/U3 34 transitions (see Fig. 5 and the 
discussion in lEllison et alj|200lh . This suggests that the stronger 
21cm component does not arise along the line-of-sight to the opti- 
cal core. While it is possible that the component arises against the 
extended radio emission that is resolved out in the VLBA image, 
this would require the absorbing clouds to be either very smooth 
and extended or, if compact, to have an extremely high 21cm op- 
tical depth. The more likely situation is that one or more compact 
CNM clouds cover part of the radio core, (but not the optical core) 



and give rise to this 21cm component. In contrast, the weaker 21cm 
component lies close to the deepest FeII-/U 122 and CII* absorption 
and fairly strong metal-line absorption is seen over most of its ve- 
locity spread. This is consistent with at least part of the latter 21cm 
component arising from gas that does lie in front of the optical 
emission. The absorbing gas thus indeed appears to have structure 
on the scale of the radio core. 

Despite the presence of this small-scale structure, it can be 
shown that if the HI is indeed in a two-phase medium as in the 
Galaxy, the WNM phase is likely to be dominant for the ab- 
sorber as a whole, on the scales probed by the radio observa- 
tions. This is because the HI in such two-phase models consists 
of randomly-distributed dense CNM "clouds" surrounded by a dif- 
fuse, widespread WNM envelope. If we assume that the optical 
QSO lies along a line of sight for which the CNM fraction is 
~ 50% (i.e. the scenario favoured by W03), the remaining half of 
the HI column density along this sightline is WNM, i.e. A^wnm ~ 
9 x 10 20 cm -2 . Th e number density in stable WNM is < 0.1 cur 3 
dWolfireetal.1 19951) ; a WNM cloud with N m ~ 9 x 10 20 cnr 2 must 
thus have a linear extent of ~ few kpc along the line-of-sight, far 
larger than the VLBA upper limit on the size of the radio core. 
The entire radio core must hence be covered by the WNM, with 
an average WNM column density of ~ 9 x 10 20 cirT 2 . Given the 
negligible 21cm optical depth of the WNM, the detected 21cm ab- 
sorption components of Fig. [TjD] must arise in the CNM. Using 
a nominal CNM spin temperature of ~ 100 K, the average CNM 
column density over the radio core (i.e. the core detected in the 
VLBA image) is Nam = 1.823 x 10 18 x 100 x (1/0.69) x 0.714, 
i.e. M:nm ~ 1.9 x 10 2C cirT 2 . The average CNM fraction across the 
radio core is hence A^cnm/C^cnm + ^wnm) ~ 0.17. Thus, even if 
the line-of-sight to the optical QSO has a fortuitously high CNM 
fraction, the large-scale HI distribution must be dominated by the 
WNM. 

We note that it is still possible to have a large fraction of the to- 
tal HI in the CNM phase if the CNM were to be distributed in small 
"clouds" with velocity separations larger than the individual veloc- 
ity widths, so that the absorption from each cloud lies below the 
GMRT 21cm detection threshold. As always for absorption stud- 
ies, such a model, in which most of the CNM is "hidden", cannot 
be ruled out from the absorption data alone. 



5 SUMMARY 

In summary, we report the GMRT detection of 21cm absorp- 
tion from the z ~ 3.39 DLA towards PKS 0201+113, with the 
21cm profile consisting of two separate components spread over 
— 1 15 km s _1 . The stronger of these components is in excellent 
agreement i n both redshif t and F WHM with the feature tentatively 
detected by iBriggs et al.1 Jl997l) with the Arecibo telescope. We 
suggest that the higher peak optical depth measured by the latter 
could be due to a combination of an incorrect flux density scale 
and low-level RFI in the single-dish spectrum, resulting in an over- 
estimate of the line depth. We find no significant evidence for opti- 
cal depth variability in the 21cm line. 

We obtain a covering factor of ~ 0.69 from a 328 MHz VLBA 
image and use this to estimate an average spin temperature of 
T s ~ [(955 ± 160) x (//0.69)] K, assuming that the average column 
density across the radio core is the same as that measured from the 
Lyman-a line. A comparison between the optical low-ionization 
metal-line and 21cm profiles finds evidence for structure in the ab- 
sorbing gas across the radio core, bringing the above assumption 



HI 21cm absorption at z ~ 3.39 9 



into question. Despite this, the GMRT observations indicate that 
most of the HI towards the radio core must be in the warm phase, 
even if a high CNM fraction is present on the narrow sightline to- 
wards the optical QSO. Finally, the complexity of the 21cm and 
CII* absorption lines and the large offset between peak optical and 
21cm redshifts suggests that the z ~ 3.39 DLA is unlikely to be of 
use in probing fundamental constant evolution. 
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